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^ ' Abstract. In this work we explore the possibility to perform "effective energy" studies in very 

, high energy collisions at the CERN Large Hadron Collider (LHC). In particular, we focus on 

^^«0 i the possibility to measure in pp collisions the average charged multiplicity as a function of the 

effective energy with the ALICE experiment, using its capability to measure the energy of the 
leading baryons with the Zero Degree Calorimeters. Analyses of this kind have been done at lower 
On ' centre- of-mass energies and have shown that, once the appropriate kinematic variables are chosen, 

, particle production is characterized by universal properties: no matter the nature of the interacting 

1^ ■ particles, the final states have identical features. Assuming that this universality picture can be 

extended to ion-ion collisions, as suggested by recent results from RHIC experiments, a novel 
approach based on the scaling hypothesis for limiting fragmentation has been used to derive the 
expected charged event multiplicity in AA interactions at LHC. This leads to scenarios where the 
^ , multiplicity is significantly lower compared to most of the predictions from the models currently 

■ used to describe high energy AA collisions. A mean charged multiplicity of about 1000 — 2000 per 

rapidity unit (at r] ^ 0) is expected for the most central Pb — Pb collisions at y/sJjN — 5.5 TeV. 
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1 Introduction 



In high-energy particle coUisions, bulk event properties 
like the average charged particle multiplicity are regarded 
as experimental observables of fundamental interest, pro- 
viding important information on the dynamics of the in- 
teraction. In particular, the average charged particle mul- 
tiplicity in multihadronic final states has so far been mea- 
sured in many different interaction systems (e+e" and 
pp{p) collisions, DIS processes, etc.) and over a wide range 
of centre-of-mass energies. Although the data show a de- 
pendence on ^/s which is characteristic of the specific ini- 
tial state under consideration, as pointed out in [UllKS] a 
universal behaviour can actually be identified if the appro- 



priate definition of the energy available for particle pro- 
duction (the "effective energy") is used. 

The aim of this work is to address the possibility to 
perform an effective energy study at energies of several 
TeV at LHC, with the ALICE experiment [l[5li]- With 
this respect ALICE (A Large Ion Collider Experiment) has 
an excellent capability, thanks to the presence of several 
detectors for the measurement of the particle multiplic- 
ity over a wide rapidity range [3[8]. Moreover, on both 
sides of the beam interaction point, the detector will be 
equipped with very forward calorimeters, the Zero Degree 
Calorimeters (ZDCs) |9], which will allow to derive the 
effective energy on an event-by-event basis by measuring 
the energy of the leading nucleons. 
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This paper is organized as follows. In section [2] we 
briefly review the main experimental results which sup- 
port the existence of a universal behaviour in particle 
production, independently of the nature of the colliding 
system. In section |3] a simplified analysis based on the 
PYTHIA Monte Carlo is presented, showing that this event 
generator is able to reproduce to a good extent the ex- 
perimental observations described in section [21 and can 
therefore be used to evaluate the feasibility of an effective 
energy study at the LHC energies. A more detailed anal- 
ysis is then described in section 31 where the capability of 
the ALICE detector to measure the energy of the leading 
baryons in pp collisions at the LHC is quantified using a 
realistic simulation of the ALICE ZDC detector response. 
Finally, in section \E\ we present a prediction for the to- 
tal charged multiplicity in Pb — Pb interactions at LHC, 
assuming that the universality features discussed herein 
hold also for ultrarelativistic heavy ion collisions. 

2 What we learned from previous experiments 

It is well known that the average charged multiplicity in 
e+e~ collisions follows a logarithmic dependence on the 
centre-of-mass energy ^/s. Figure [H shows a compilation 
of data from e+e" experiments [TUllTTIIT^lfnifTIllTSllTClfTTl 
fT8 l fT9 l l20 l [2T] over a wide range of centre-of-mass energies 
(full black symbols) , together with the result of a logarith- 
mic fit to the measurements, indicated by the dashed line. 
As mentioned in section [H the average charged multiplic- 
ity is characterized by a significantly different dependence 
on y/s if other initial states are considered. In particular, 
in the case of pp{p) collisions the charged particle multi- 
plicity at a fixed centre-of-mass energy is observed to be 
systematically lower than what can be inferred from the 
e+e" data at the same ^/s (see again Fig.[Tl open symbols, 
data from [ ^[^M^ill^M^ 1 ■ As pointed out in Jj, this 
behaviour can be understood after considering that, while 
in e+e~ collisions the energy available for particle produc- 
tion coincides with the full centre-of-mass energy (once 
the effects from the initial state radiation are removed), 
in pp{p) collisions this energy is reduced with respect to 
^/s due to a basic feature of the hadronic interactions, the 
"leading effect". 

The leading effect, which is related to the quantum 
number flow between the initial and the flnal state, im- 
plies that in pp{p) collisions there is a high probability 
to emit a forward baryon with large longitudinal momen- 
tum along the direction of the incident beams. Since the 
leading baryon carries away a fraction of the incident en- 
ergy, the energy available for particle production is re- 
duced with respect to the total centre-of-mass energy. If 
the leading effect is taken into account in the definition 
of the effective energy available for particle production, 
Eeff, common universal features can be observed in the 
two different interaction systems, e"'"e~ and pp{p). 

In particular, as done in [3], the effective energy can 
be estimated on an event-by-event basis by measuring 
the energy Eieading of the leading baryon in each event 
hemisphere (the hemispheres being defined with respect to 
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Fig. 1. Average charged particle multiplicity measured in 
e^e~ (full black symbols) and pp{p) (open symbols) collisions. 
A logarithmic fit to the e'^e~ data is also shown as the dashed 
line. The full grey symbols indicate the pp{p) data after the 
centre-of-mass energies have been scaled by a factor | (see 
text). 



a plane transverse to the direction of the incident beams 
in the centre-of-mass system). The effective energy per 
hemisphere (£'e//)i is then given by: 

iEeff)i = - {Eleading)i fori=l,2. (1) 

Notice that, by measuring the energy of the leading 
particles event by event, a wide range of effective energies 
can be covered at a fixed centre-of-mass energy. Rely- 
ing on the experimental observation of the independence 
of the two event hemispheres with respect to the leading 
effect 3 , the total effective energy E^f / available for par- 
ticle production in the whole event can be derived from 
the measurement of only one leading particle per event, 
using the relation: 

^eff — '2{Eeff)i ~ — 2{Eieading)i- (2) 

Correspondingly, the average total charged multiplicity in 
the event is given by : 

(rich) = 2 {nch)i . (3) 

Alternatively, in those events where a leading particle is 
measured in both hemispheres, the total effective energy 
can also be derived from the relation: 

Eeff = v/s [1 - {xf)i][1 - {xf)2], (4) 
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where {xF)i are the measured Feynman-x of the two lead- 
ing nucleons, (a;F)i = 2{pL)i/^/s. 

The two methods of defining the effective energy in the 
event have been proved to be equivalent, as shown in [1]; 
it should be noticed, however, that the first one has the 
advantage of making optimal use of the data in terms of 
the collected event statistics, obviously owing to the much 
larger acceptance and efficiency of single with respect to 
double leading particle detection. 

As mentioned before, after taking into account the 
leading effect and using the effective energy instead of 
the full centre-of-mass energy, a very good agreement be- 
tween the measured charged particle multiplicity in e"^e^ 
and pp collisions is obtained, as shown in Fig. [2] The pp 
data in this figure refer to a "minimum bias" event sample 
from which elastic and diffractive processes were excluded, 
with leading protons in the range 0.3 ^ xp 0.8 (see later 
on sections [3] and [4|) . 

A good agreement between the inclusive pp(j)) results 
presented in Fig. [1] and e+e~ results is also obtained if 
the pp{p) centre-of-mass energies are scaled by a fac- 
tor i (full grey symbols). This is a direct consequence of 
the effective energy being the relevant variable in particle 
production, and of the assumption, supported by several 
measurements [1] , that the Feynman-x distribution of the 
leading particles has a mean value equal to 0.5, result- 
ing in an average effective energy which is about half the 
full centre-of-mass energy. 
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Fig. 2. Average charged multiplicity (rich) as a function of 
the effective energy -Ee// (here indicated as \/(^£rt)^)j ^ mea- 
sured in minimum bias pp collisions collected by the SFM 
experiment at the CERN ISR (full circles). The data from 
e"'"e~ experiments are also shown (open circles and triangles) 
in terms of -^s. A fit to ISR-pp data is superimposed. The plot 
is taken from [2]. 



3 Pythia simulations 

The PYTHIA6.2f4 [27j Monte Carlo has been used to sim- 
ulate pp collisions at a centre-of-mass energy of f 4 TeV. 
An effective energy study carried on for pp collisions re- 
quires a precise definition of the type of events to be 
used. Following the choice of the original ISR-pp exper- 
iment only non diffractive minimum bias events 
have been considered herein, ft should be noted that most 
of the other pp{p) experiments [22l[23l l24''25,26 quoted 
in section [2] included double-diffractive processes in their 
definition of minimum bias. With PYTHIA the double- 
diffraction contribution is at the level of 15% with respect 
to the non diffractive one at ^/s = 14 TeV (using default 
generation parameters). A brief description of the ALICE 
pp minimum bias trigger capabilities for both diffractive 
and non diffractive processes is given in subsection 14. II 

To optimize the agreement with respect to the depen- 
dence observed in the data for the charged particle mul- 
tiplicity as a function of the effective energy, some of the 
PYTHIA parameters regulating the treatment of the mul- 
tiple parton interactions (MI) have been tuned. This MI 
model has been introduced and widely discussed in the 
literature P51I^M5D1[5T] in order to describe the basic fea- 
tures and, in particular, the charged particle multiplic- 
ity of high energy minimum bias pp{p) final states up to 
^/s ~ I TeV. It is still under investigation but represents 
however one possible tool to perform extrapolations in the 
multi-TeV domain at LHC. The set of values used in the 
simulation is reported in table [TJ 



Parameter 


Our tuning 


Default 


MSTP(51) 
MSTP(81) 
MSTP(82) 
PARP(82) 
PARP(89) 
PARP(90) 


7 (CTEQ5L) 
1 
2 

2.15 GeV 
1000 GeV 
0.215 


7 (CTEQ5L) 
1 
4 

1.80 GeV 
1000 GeV 
0.160 



Table 1. Values of the PYTHIA parameters tuned for this 
work to generate minimum bias (non diffractive) pp events at 
various centre-of-mass energies. 



Notice that the only three parameters tuned corre- 
spond, respectively, to the selection of the MI mecha- 
nism features in terms of hadronic matter distribution 
(MSTP(82)), the regularization scale of the transverse mo- 
mentum spectrum of the process (PARP(82)) which con- 
trols the MI rate and the exponent of the -yi-power law 
dependence of the MI mechanism (PARP(90)) 

Figure [3] shows the mean total charged multiplicity for 
generated pp collisions over a wide range of centre-of- 
mass energies, with ^/s rescaled by a factor ^ to take into 
account, on average, the energy carried away by the lead- 
ing nucleons, as done in section [2] To derive the average 
charged multiplicity {rich) all primary charged particles in 
the event (including the decay products of short-lived res- 
onances) with no transverse momentum cutoff and in the 
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full angular acceptance have been considered. Although 
{rich) is derived from the Monte Carlo information at the 
event generator level, notice that the ALICE detector is 
expected to measure the multiplicity with a good accu- 
racy, high efficiency and over a wide acceptance both in 
case of pp and Ph — Ph collisions (as reported in chapter 
6.1.5 of the ALICE Physics Performance Report vol. II 

m)- 

Using the tuned PYTHIA Monte Carlo, a very good 
agreement with the fit that describes the y/s dependence 
of the multiplicity for e+e" collisions (see Fig. [T]) is ob- 
served in Fig. [3] over the whole range of energies. Using 
instead the default values, PYTHIA disagrees by ~ 10% 
with the e+e~ fit for > 1 TeV. 
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Fig. 3. Mean charged multiplicity {rich) in pp collisions from 
PYTHIA Monte Carlo, at ^ = 70 to 14000 GeV. The energy 
on the horizontal axis is scaled by a factor ^ in order to take 
into account an average leading effect in the pp final states. 
The e"'"e" fit of Fig. [T]is also plotted. 



A satisfactory agreement is also observed when study- 
ing the multiplicity as a function of the effective energy in 
each single event. Following the approach described in sec- 
tion [1] for each event the effective energy in a single hemi- 
sphere was derived from the energy of the corresponding 
leading baryon (defined as the highest longitudinal mo- 
mentum baryon in the hemisphere) and then multiplied by 
2, according to relation ([2]). The total charged multiplicity 
was also derived on a single hemisphere basis, scaled by 
the same factor of 2. As can be seen from Fig. 01 also in 
this case the dependence of the total charged multiplicity 



on the effective energy follows closeljQ the fit to the e+e 
data. 

Notice that in the above analysis the Feynman-x of 
the leading baryons was restricted to be in the range: 
0.3 < < 0.8. The lower cutoff, a;^™ = 0.3, is mainly 
motivated by the fact that in a realistic analysis the iden- 
tity of the leading particle will not be accessible. Some care 
is therefore required in defining an appropriate Feynman- 
X acceptance, in order to ensure that the energy measured 
in the forward region is indeed in most cases due to a 
leading baryon, and that the contamination from leading 
particles which are not baryons is kept below a reasonably 
small level. In Fig. [5] the xp distribution of the leading 
baryon, defined as the proton/neutron with the highest 
longitudinal momentum (i.e. with the highest Feynman- 
X = 2pL / y/s) among all baryons produced in a given 
event hemisphere, is shown as the dark-shaded histogram, 
for 10^ pp collisions generated at ^/s = 14 TeV. The xp 
distribution for the cases where this leading nucleon is 
also the leading particle among all particles (baryons or 
non baryons) of the same hemisphere is superimposed as 
the light-shaded histogram. For Feynman-x greater than 
0.3, it can be seen that the leading nucleon is indeed the 
leading particle in the event hemisphere. Moreover, a re- 
quest of a minimum xp allows to reduce the contamination 
from leading particles that are not baryons. At xp = 0.3 
this contamination is ~ 30% and rapidly decreases for 
xp > 0.3. 

For what concerns the constraint on the maximum xp, 
the requirement a;™°^ = 0.8 was ispired by the original ef- 
fective energy data analysis [Il|2l[3], where the same kind 
of cutoff on leading protons was imposed. This xp < 0.8 
condition, which fairly matches the acceptance of the AL- 
ICE ZDC calorimeters for leading baryons (see later on 
section 131), would also imply the rejection of the protor0 
single-diffraction peak (not considered in this study). 

The PYTHIA generator reproduces to a good extent 
also other features related to the leading effect which have 
been observed in the data. As an example. Fig. I6l shows, 
for pp collisions at y/s ~ 14 TeV, the ratio between the 
total charged multiplicity (rich), measured as a function of 
the total effective energy Eeff given by equation and 
the average charged multiplicity in a single hemisphere, 
measured as a function of the effective energy in a sin- 
gle hemisphere, scaled by a factor 2, Eeff = 2{Eef f)i,2- 
The ratio is consistent with the relation (rich) = 2 {uchji 2 
in agreement with the hypothesis of the two event hemi- 
spheres being independent with respect to the leading ef- 
fect, and of the equivalence of the two effective energy 
Eeff definitions presented in section [l] equations Q and 



^ The values deviate systematically from the e^e~ data fit 
as Ecf f increases whereas they always start well on the fit line. 
This effect is likely due to the fact that the multiple interaction 
parametrization of PYTHIA generator slightly depends also on 
the centre-of-mass energy. 

^ For neutrons single-diffraction processes with charge ex- 
change are highly unlike. 
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Fig. 4. Mean charged multiplicty {rich) in PP PYTHIA events 
at ^/s = 70 to 14000 GeV with respect to effective energy. The 
e'^e~ fit is also plotted. Here the single hemisphere variables 
are considered and scaled by a factor of 2 (see text). 
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Leading Particle = Baryon 




Fig. 5. The xf distribution of the most energetic baryon 
(dark-shaded histogram) in each hemisphere for pp collisions 
at -^s = 14 TeV. The xf distribution of this baryon in the 
case it is also the most energetic particle in the hemisphere is 
shown as the light-shaded histogram. 



4 Tagging of leading nucleons with the 
ALICE ZDC detector in pp collisions 

Although the programme of the ALICE experiment has 
its main focus on heavy ion physics, the detector will also 
have excellent physics capabilities in the case of pp col- 
lisions. In particular, this applies to the analysis we are 
presenting in this paper. 



4.1 ALICE performance for the selection and 
measurement of pp events 

A detailed description of the ALICE trigger system can 
be found elsewhere [5115^. Here we will limit ourselves to 
a brief report on both the minimum bias trigger selection 
and the charged particle multiplicity measurement [8 that 
ALICE will be able to perform in pp collisions at LHC. 

The minimum bias trigger is realized coupling different 
signals coming from the VO [33] and the SPD [34] detec- 
tors. The general purpose of this minimum bias trigger is 
to select events from pp collisions with an cfliciency as high 
as possible, and a bias as low as possible, and to reject 
events due to beam-gas or beam-halo. The VO is made 
up of two different arrays of scintillators placed along the 
beam pipe in the forward/backward directions. It provides 
two kinds of signals: (i) VZERO.OR that requires at least 
one hit in one counter on one side and (ii) VZERO.AND 
that requires at least one hit in one counter on both sides. 



The SPD (Silicon Pixel Detector) is made up of two coax- 
ial layers of silicon sensors which also provide a trigger 
signal (GLOB.FO) consisting in a fast-OR signal of all 
its 1200 chips (400 in the inner layer, 800 in the exter- 
nal one). Different combinations of these signals, coupled 
to a non background trigger (notBG, defined by requir- 
ing no signals in either of the VO counters within appro- 
priate time windows corresponding to beam-background 
processes), set different trigger configurations. The three 
configurations are listed in table [21 while the correspond- 
ing efhciencies are reported in table [3j 



MB 1 


GLOB.FO 


or 


VZERO.OR 


and 


notBG 


MB 2 


GLOB.FO 


and 


VZERO.OR 


and 


notBG 


MB 3 


GLOB.FO 


and 


VZERO.AND 


and 


notBG 



Table 2. Signal configurations for pp minimum bias trigger. 



The multiplicity can be measured counting the "track- 
lets" (extracted by the association of the signals in both 
the two layers of the SPD) with the SPD in the central 
region and with the FMD (Forward Multiplicity Detec- 
tor) [33] in the forward region. The pseudorapidity range 
covered by the ALICE detector is about 8 7y-units. In this 
range an accurancy on the multiplicity better than 10% 
with 0.1 77-units bin width can be achieved (as estimated 
in single central HIJING [35] events). 
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Fig. 6. Average charged multiplicity ratio (rich) / {wch)i as a 
function of the total effective energy E^f / (black symbols) and 
average charged multiplicity ratio {rich) 2 I (^ch)i as a function 
of the effective energy in a single hemisphere scaled by a factor 
2, -Be// = 2(£'e//)i,2 (grey symbols). 



ment consists of two distinct calorimeters, one for leading 
neutrons (ZDCN), placed at zero degrees relative to the 
LHC axis, and one for leading protons (ZDCP), placed ex- 
ternally to the outgoing beam pipe on the side where pos- 
itively charged particles are deflected by the beam optics. 
For charged particles, the Feynman-x range experimen- 
tally accessible with the ZDCs will be G [0.30 - 0.64], 
as the beam optics selects energies between 2.2 TeV and 
4.5 TeV. For neutral particles such constraints are clearly 
not present. A more detailed description of the calorime- 
ters is given in [3]. 



Neutron calorimeter 




Vocuum pipe aperture 068( 

Vacuum pipe aperture 
optional shape 



Beam Pipe 044/40 (ring 2) 
Proton calorimeter 



Process (ctmtcV [mb]) MBf (%) MB2 (%) MB3 (%) 
Non Diff. (55.22) 99.9 99.1 96.9 
Single-Diff. (14.30) 73.8 59.5 38.4 
Double-Diff. (9.78) 87.8 68.7 45.6 



Fig. 7. Transverse section of the LHC beam line at a distance 
of 115 m from the interaction vertex. The location of the two 
ZDC calorimeters for protons and neutrons is also shown. 



Table 3. Minimum bias trigger efficiency for different (diffrac- 
tive and non diffractive) pp processes. 



In the following we will assume that the uncertainty 
on the multiplicity measurement is negligible with respect 
to that on the effective energy measurement. 

4.2 ZDC performance for leading nucleons 

As mentioned in section [T] the ALICE experiment will be 
able to measure the energy of the leading particles needed 
for this study thanks to the presence of the Zero Degree 
Calorimeters (ZDCs) [9] . The aim of this section is to show 
that this system, although designed for the detection of 
nearly monochromatic {Ef^ « 2.75 TeV) leading nucleons 
in ion-ion collisions at ^/sJin = 5.5 TeV, will be able to 
measure on an event -by-event basis leading particles in pp 
collisions over a wide range of energies and with a good 
accuracy. 

A scheme of the ZDC apparatus is shown in Fig. \7\ it 
consists of two identical elements, placed on both sides of 
the interaction vertex at a distance of 115 m. Each ele- 



The performance of the ALICE ZDC calorimeters were 
first checked on single highly energetic protons and neu- 
trons, with generated longitudinal momenta in the range 
between 1 TeV and 6 TeV. As mentioned before, due to 
the beam optics the ZDC acceptance is restricted to pro- 
tons with energies in the range 2.2 < E < 4.5 TeV, while 
in the case of neutrons the whole range of energy is acces- 
sible. The GEANT '36] package and the ALICE simula- 
tion and reconstruction software, AliRoot [37], were used 
for the simulation of the detector response, and to per- 
form realistic digitization and reconstruction in the ZDC 
calorimeters. 

In Fig.[5]the energy reconstructed in the ZDCs for pro- 
tons and neutrons as a function of the generated energy 
is shown. While a good linearity is observed in the case of 
neutrons over the whole range of energies, in the case of 
protons the reconstructed energy is underestimated with 
respect to the true energy above ~ 3 TeV. This is because 
with increasing energy, the impact point of these particles 
falls closer and closer to the edge of the detector due to the 
beam optics, and the hadron cascade is no longer fully con- 
tained in the calorimeter. The energy resolution, shown in 
Fig. [HI ranges between 10% and 20% both in the neutron 
and proton cases in the selected energy range (from 2.2 to 
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and neutron (black symbols) ZDC calorimeters, as a function 
of the proton/neutron generated energy. 
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4.5 TeV). The absolute error on the energy measurement 
is about 500 GeV. For protons the ZDCP energy resolu- 
tion could likely be improved using the information of the 
impact point of the track on the calorimeter front face. 

After checking that the ZDCs can measure the energy 
of forward-going nucleons over an extended xp range, 



their performance has been tested on fully simulated pp 
collisions. Starting from 14000 generated events, in about 
75% of the events a signal is recorded in either the proton 
(ZDCP) or the neutron (ZDCN) calorimeter, in the for- 
ward or backward regions. In the majority of cases the pro- 
ton calorimeter is hit by a single particle, while in the case 
of the neutron calorimeter the multiplicity is higher, due to 
additional photons coming from the decay of high-energy 
forward tt'^ mesons produced in the interaction (in the case 
of the proton calorimeter, the corresponding contamina- 
tion which would be expected from charged pions is to a 
large extent removed by the beam optics). This is demon- 
strated in Fig.sHnia) and [TOlb), where the multiplicity of 
particles detected in each of the ZDC calorimetertl3 and 
their particle species composition are respectively shown. 
Moreover, looking at the fraction of deposited energy of 
the most energetic particle hitting the ZDCs (Fig fTTj) . it 
can be seen that in nearly all cases more than 95% of the 
total energy released in the ZDCP is due to a single par- 
ticle, while in the ZDCN this happens only in about 50% 
of the cases. This is again mostly due to the background 
coming from photons from 7r° decays mentioned above. 

The presence of additional particles impinging on the 
neutron calorimeter causes a substantial distorsion when 
the total energy measured in the calorimeter is used to 
estimate the energy of the leading particle entering the 
ZDCN, as can be seen in Fig[T2lb, where the difference be- 
tween the total energy measured in the ZDCN calorimeter 
and the generated leading particle energy is shown. 

In the case of the proton calorimeter, where this type 
of background is much lower, the distribution is nearly 
gaussian (see Fig[T2la), although the energy is, on average, 
slightly underestimated due to hadronic showers not being 
fully contained for E > 3 TeV. Therefore, in terms of 
the effective energy reconstruction, the proton calorimeter 
ZDCP is expected to provide a cleaner measurement with 
respect to the ZDCN. On the other hand, the neutron 
calorimeter covers a range in Feynman-x which is more 
extended than in the case of the proton calorimeter, as 
discussed before in this section. 

The reconstructed x^p"^" = IEzdc I distribution, 
where Ezdc is the total energy in cither the ZDCP or 
the ZDCN, is shown in Fig [HI 

In those cases where a signal is detected in both ca- 
lorimeters, the highest among the two energies was taken 
as an estimate of the leading particle energy. The empty 
and the dark-shaded histograms indicate respectively the 
Feynman-x distribution of the reconstructed leading par- 
ticles and the contamination from cases where the highest 
energy deposit in the calorimeter was not due to a leading 
baryon. This latter component, which can be regarded as 
a source of contamination in the measurement of the ef- 
fective energy, is about 30% in the case of the ZDCP and 
40% in the case of the ZDCN over the Feynman-x range 
used in the analysis (0.3 < x^/'^° < 0.8). This contami- 
nation, however, is expected not to affect significantly the 
quality of the measurement of the multiplicity as a fuction 



In order to reject very low-energy background, only parti- 
cles with Feynman-x xf > 0.1 were considered. 
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Fig. 12. Difference between the reconstructed and the gener- 
ated energy in the ZDCP (a) and ZDCN (b). In the case of 
the ZDCN, the clear distortion on the positive side of the dis- 
tribution in (b) is due to the presence of high energy 7's from 
tt" decays, which cause an overestimate of the reconstructed 
energy. 



of the effective energy. This is shown in Fig. [TH where the 
generated average charged event multiplicity as a fuction 
of the reconstructed effective energy E^jy = y/s — 2EzDC 
(black symbols), and of the effective energy calculated at 
the generator level (grey symbols) are compared. In both 
cases equation ([2]) was used. The width of the horizontal 
binning is determined by the energy resolution of the ZDC 
calorimeters. 

The results are consistent, indicating that the limita- 
tions related to the reconstruction of the effective energy 
in the ZDCs and the background from other particles in 
the forward region do not introduce a significant bias in 
the measurement. We therefore expect that the ALICE 
experiment will have a good capability to perform an ef- 
fective energy study in collisions at LHC. 



5 Nucleus-nucleus collisions 

As described in the previous sections, hadron-hadron col- 
lisions produce a mean charged multiplicity in good agree- 
ment with the e~ collisions at the same effective energy. 
Another question is still open. What happens in nucleus- 
nucleus? The aim of this section is to investigate the pos- 
sible scenario which we expect to be observed at LHC. 

The PHOBOS experiment at RHIC has recently shown 
(Fig. [T5)) how the total charged multiplicity in AA colli- 
sions, scaled by the number of participant nucleons, varies 
with energy |38J. In particular, the AA data indi- 

cate that the scaled total charged multiplicity does not 
depend, to a good approximation, on the centrality of the 
collision and that its value vs. yT/vw, for 20 < y^sjvw < 
200 GeV/nucleon pair, is in good agreement with the one 
measured in e^e~ collisions vs. ^/s. These findings favour 
the hypotesis that in A A collisions, at very high energy. 
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Eeff would coincide with y^JJvTV for each nucleon-pair 
system. 

Assuming that this will hold also at LHC energies, it 
is then possible, on the basis of an extrapolation from the 
e^e~ fit, to derive a prediction for the mean total charged 
multiplicity per nucleon pair in nucleus-nucleus collisions 
at yfsNN = 5.5 TeV. The predicted mean total charged 
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Fig. 15. (a) Average charged multiplicity for pp, pp, e e~ 
and central Au — Au events as a function of the centre-of- 
mass energy per binary collision. The Au — Au {Pb — Pb for 
NA49 experiment) data are scaled by a factor {Npart/'2). The 
dotted line is a perturbative QCD fit to the e^e'^ data. The 
diamonds are the pp{p) data plotted at energy values scaled by 
a factor 1/2. (b) The data in (a) divided by the e"'"e~ fit, to 
allow direct comparison of difi'erent data at the same energy. 
Both plots (a) and (b) are taken from [38j . 



multiplicity, scaled by the number of participants (Npart), 
would in fact be: 



{nch) 



72 ± 2, 



part 



(5) 



resulting then in a total charged multiplicity in very cen- 
tral Pb — Pb collisions at the LHC (using the maximum 
possible value Npart = 416, which is the total number of 
nucleons available in Pb — Pb collisions): 



{rich) = 15000 ±400. 



(6) 



Once the total multiplicity is fixed, a prediction for 
the charged multiplicity at midrapidity (which is relevant 
for the measurement of a lot of observables in Pb — Pb 
collisions) can be derived in the so-called limiting frag- 
mentation hypothesis. As shown in Fig. [TBI several mea- 
surements from pp{p) [39ll40] at different centre-of-mass 
energies support the fact that the multiplicity at rapidities 
close to the beam rapidity (the so-called fragmentation re- 
gion) has a shape which is independent of the energy in 
the collision. This effect [H] has been observed at very 
different centre-of-mass energies also in DIS events. It 
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depends only on the nature of the incident systems and, 
in the case of nuclei, on the centrality of the collision. 
The same behaviour is indeed preserved in AA collisions 
at RHIC, as shown in Fig. [17] by the measurements of 
the rapidity distribution for different centralities from the 
PHOBOS collaboration. 



Therefore, in these hypotheses it is possible to fix the 
shape of the pseudorapidity distribution [M1|401I41II42II43[ 
44 , 45 and derive the mean charged multiplicity at midra- 
pidity in central Pb — Pb collisions at LHC, as it is shown 
in Fig. [m 




^■ybeam 

Fig. 16. Limiting fragmentation in pp and pp collisions at 
different centre-of-mass energies. 
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Fig. 18. Limiting fragmentation at LHC based on an extrap- 
olation from RHIC results in central (0-6%) ion-ion collisions. 
The dashed line indicates the value of the charged multiplicity 
per rapidity unit (scaled by the number of participants) which 
results in a total multiplicity (given by the integral under the 
curve) equal to the prediction from the e'^e~ fit. The rapidity 
of the beam at LHC is yteam ~ 8.6. 
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Fig. 17. PHOBOS (RHIC) results for Au - Au collisions at 
y/SNN = (19.6, 200) GeV. In this figure, taken from [S^, the 



-dNch/drj' , where rj' = rj — ybeam, are presented 



results of 



for different centrality bins. 



As shown in [39ll40ll4T] . the behaviour of the pseudo- 
rapidity distribution has some universal features and this 
allows to fix the shape of the pseudorapidity distribution a 
priori, independently of the centre-of-mass energy of the 
collision. Anyway such universal behaviour doesn't allow 
to establish the height at which the distribution is cut. 
This information can be found requiring that the integral 
under the curve should give the total charged multiplicity 
as predicted from the extrapolation logarithmic fit to the 
e~^e~ data. 



The result is: 
2 



part 



duch I 



5.3 ± 0.5, 



(7) 



corresponding to about 1100 ± 100 charged tracks per 
rapidity unit. 

There are a lot of predictions for this observable based 
on as many models. These models, that are compatible 
with the actual experimental results or are tuned on them, 
produce a very wide spread of predictions. In particular, 
a recent direct extrapolation from RHIC data seems to 
indicate that the charged multiplicity at LHC, in heavy 
ion central collisions at midrapidity, could be very small 
{{drich/dri) ~ 1200) [46], in good agreement with our pre- 
diction. 

Figure [TH] taken from [ITj is a review of the results of 
the main models. It can be noted that there are predictions 
of over 6k charged tracks per rapidity unit and others of 
less than 3k. 

Our prediction is significantly low when compared with 
the predictions of widely used Monte Carlo generators 
(such as HIJING, for example). It should be noted how- 
ever that other effects, in particular the jet quenching |48j . 
could sizeably increase the multiplicity in the central ra- 
pidity region. In Fig. [19] only the point relative to HIJING 
takes somehow into account such an effect (averaging on 
the two possible scenarios with and without quenching) 
[?7] . An increase by a factor ~ 2 is actually predicted 



11 



dN/dy at y = for charged particles at LHC, 1 /20 centra lity 



10 12 14 16 



Fig. 19. Mean charged multiplicity per rapidity unit at r; = 
from different models, at ^/sNN = 5.5 TeV. For those mod- 
els for which the extrapolation to LHC energies has not been 
performed, the corresponding value is set to zero [47] . 



by HIJING when its defaullQ jet quenching algorithm is 
switched on. Even considering this additional effect, if the 
universality features hold at LHC energies, a relatively 
low value for the charged multiplicity at midrapidity is 
expected to be observed. In that case, while the physics 
reach of ALICE will not be endangered, as illustrated in [5] 
versus many representative physics observables, its back- 
ground conditions would correspondingly improve. 



6 Conclusions 

In this work we have estimated the average charged par- 
ticle multiplicity at midrapidity of the events that will be 
produced in Pb — Pb interactions at LHC, in an extremely 
high energy density domain. Our result, derived by means 
of detailed Monte Carlo simulations, was inspired by pre- 
vious experimental findings and based on the sound hy- 
pothesis that the universality features of multihadron final 
states will still apply at LHC, in both pp and A A interac- 
tions. 

The measurement of the leading baryons which will 
be needed to deduce the effective energy, e.g. the essen- 
tial parameter of this analysis, has been proven to be fea- 
sible using the Zero Degree Calorimeters (ZDCs) of the 
ALICE detector. We have shown, in particular, the capa- 
bility of ALICE to perform an effective energy study in 
the range 3—10 TeV for pp collisions at ^/s = 14 TeV. 



In ALICE different variants of this algorithm are presently 
being investigated. 



This will be possible taking advantage of the acceptance 
and energy resolution (~ 500 GeV) of the ZDC calorime- 
ters to measure protons and neutrons with Feynman-x 
0.3 <XF < 0.8. 

The rather low multiplicity prediction obtained herein 
for large centrality Pb — Pb collisions at LHC, namely 
{duch/dri) = 1000 - 2000 at - and ,JIWn = 5.5 TeV, 
adds to the so far uncertain mosaic of the predictions of 
such a very fundamental quantity for the description of 
AA phenomena. 
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